Self-assembly of amphiphilic peanut-shaped nanoparticles 
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We use computer simulation to investigate the self-assembly of Janus-like am- 
phiphilic peanut- shaped nanoparticles, finding phases of clusters, bilayers and mi- 
celles in accord with ideas of packing familiar from the study of molecular surfactants. 
However, packing arguments do not explain the hierarchical self-assembly dynamics 
that we observe, nor the coexistence of bilayers and faceted polyhedra. This coexis- 
tence suggests that experimental realizations of our model can achieve multipotent 
assembly of either of two competing ordered structures. 

Components are said to 'self-assemble' when they organize to form stable patterns or 
aggregates without external direction. Self-assembly is driven by interactions as differ- 
ent as weak covalent bonds and capillary forces, involves components ranging in size from 
Angstroms to centimeters, and occurs both in inorganic settings and in living organisms p]- 
[9]. Mimicry of the self-assembly seen in the natural world promises the development of new, 
functional materials patterned on the nanometer scale [10l HI]. In pursuit of this goal, we 
take inspiration from the self-assembly of molecular surfactants to investigate using com- 
puter simulation the behavior of a simple model of their colloidal counterparts. Molecular 
surfactants, comprising chemically linked hydrophobic and hydrophilic groups, are of central 
importance in biology and industry, able to form a plethora of phases in water or mixtures 
of water and oily liquids. These phases include micelles [12], bilayers and vesicles, as well as 
numerous bicontinous phases [13] that serve as internal cellular packaging [14] and are the 
bane of many a plumber. 

Here we ask: What might self-assemble in aqueous solution from amphiphilic, peanut- 
shaped colloidal nanoparticles? Such particles can now be prepared in large quantities, 
starting from spherical, crosslinked polystyrene 'seed' particles of radius ^ 50 nm. Mixing 
these seeds with styrene initiates monomer-polymer phase separation and creates particles 
with a peanut-like shape characterized by two fused spherical lobes of controllable relative 
size. Additional treatment of the seed particle renders peanuts amphiphilic, with one lobe 
hydrophobic and the other hydrophilic [T5HT7] : the resulting body can be regarded as a 
generalized Janus particle [18J . (We note that colloidal silica dumbbells can be synthesized by 
other routes [19], and that micrometer-scale peanuts show potential as Pickering stabilizers of 
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oil-in- water emulsions [20J.) In an attempt to answer our posed question we have constructed 
a model of interacting peanuts whose minimal character is motivated by the insight into 
self-assembly afforded by similarly simple model systems |2lT - [39| . Our model can be evolved 
with computational efficiency sufficient to allow observation of collective, thermally-driven 
dynamics on timescales of seconds. Its construction rests upon two assumptions: first, 
that peanuts in aqueous solution experience a thermodynamic driving force that causes 
hydrophobic lobes to attract each other; and second, that functionalization of peanuts' 
hydrophilic lobes renders them chemically passive. We assume solution conditions to be 
such that electrostatic interactions between peanuts mediate only short-ranged repulsions. 

Model geometry is shown in the Appendix (Fig. [5]). Each peanut consists of a spherical 
hydrophobic lobe of radius R fused with a spherical hydrophilic lobe of radius R\. The 
centers of the two lobes are separated by a distance 5, which we quantify in dimensionless 
form via the parameter e = S/(Rq + R\). We fixed R throughout (we consider R to 
be approximately 50 nm), and investigated the behavior of the model as we varied e and 
Ri/Rq. We required \Ri — R \/(Ri + R ) < e < 1; the lower inequality stipulates that one 
lobe may not be completely buried in the other, while the upper inequality requires lobes 
to be in contact. We define the orientation vector Si of peanut i to be the unit vector 
whose origin is the center of the hydrophobic lobe and which points diametrically away from 
the center of the hydrophilic lobe. Neighboring peanuts interact via hard-core excluded 
volume interactions - nothing may approach closer than R (resp. i?i) to the center of 
each hydrophobic (resp. hydrophilic) lobe - and via a short-ranged, pairwise interaction 
modeling the attraction between solvated nanoscale hydrophobes [40] (we do not consider 
solvent explicitly in our model). We assume the hydrophobic interaction to be attenuated on 
a scale of approximately 5 nm, and impose an attractive interaction between the hydrophobic 
lobes of peanuts i and j of the form 

Un = e b 6(r c - r l3 )A(e i3 )A{e 3l ) (C a (r ZJ ) - £ a (f c )) . (1) 

Here is a binding energy; x = x/R — lisa shifted and scaled distance; ry is the distance 
between the centers of the hydrophobic lobes of i and j, respectively Q and Cj] r c = 2.5i?o 
is a cutoff length; and C a (x) = A(x~ 2oL — x~ a ) is a generalized Lennard- Jones function. We 
take a = 15 to ensure an attraction of sufficiently short range. In simulations we varied 
the attractive binding energy between limits of 4/cbT and 8/cbT: corresponding forms of 
the radial component of the hydrophobic interaction potential are plotted in Fig. [5](a). In 
experiment, variation of the strength of the hydrophobic driving force may be achieved by 
variation of temperature or solvent composition. 

The factors A in Equation ([T]) parameterize an angular modulation of the attractive 
interaction, with being the angle between Si and the vector pointing from Q to Cj. This 
modulation enforces the tendency of two hydrophobes of girth exceeding 1 nm to maximize 
their surface-to-surface contact [40]. Its form is given and derived in the Appendix. 

We performed simulations of collections of peanuts of given geometry, defined by their 
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values of Ri/Ro and e, for a range of values of the hydrophobic attraction strength e\>. We 
used 1000 peanuts in a periodically-replicated cubic simulation box of side 64i? , correspond- 
ing to fixed mole fraction of colloid (this choice reflects the low concentrations that we intend 
to use in experiments; in exploratory simulations we verified that the phases seen here also 
form at several other choices of (low) concentration). The bulkiest (Ri/Rq = 2, e = 1) and 
most compact (Ri/Rq = 0.2, e = 0.8) peanuts we considered occupied volume fractions of 
about 14% and 1.5%, respectively (see Appendix). Starting from configurations in which 
peanuts were randomly mixed and oriented, we evolved each system according to the ver- 
sion of the virtual-move Monte Carlo algorithm [41] described in the Appendix of Ref. [42] . 
This algorithm is designed to mimic an overdamped dynamics by making collective moves of 
particles according to the potential energy gradients, or forces, they experience; such collec- 
tive moves are neglected by standard single-particle Monte Carlo algorithms. In brief, one 
particle is selected and subjected to a translation or a rotation. If changes in pairwise poten- 
tial energies between that particle and its neighbors are favorable then the chosen particle 
moves independently; if not, neighbors are recruited iteratively and experience a collective 
displacement or translation. The acceptance rate for each move is chosen to preserve de- 
tailed balance and to reflect, in an approximate fashion, the (short-ranged) hydrodynamic 
drag suffered by the collective body. We scaled collective translation acceptance rates by 
the reciprocal of an approximate hydrodynamic radius of the moving body, appropriate for 
its instantaneous direction of motion [41] (identical for forward and reverse moves, as re- 
quired to satisfy detailed balance), and scaled acceptance rates for rotations by the cube of 
a similar factor, accounting for an aggregate's instantaneous axis of rotation (identical for 
forward and reverse rotations). While these damping factors are approximate, control of 
aggregate diffusion constant scalings - and particularly enforcement of the notion that an 
aspherical body does not translate equally rapidly in all directions - provides one advantage 
of this method over conventional integration of overdamped equations of motion. We drew 
translation magnitudes from a uniform distribution with maximum 0.3i?o [55j, and drew 
rotation angles from a uniform distribution with maximum 13.6° (rotations were performed 
about a randomly-chosen axis through the midpoint of the line joining the centers of hy- 
drophobic and hydrophilic lobes). This ability to make large trial translations and rotations 
of individual particles in the face of attractions and repulsions that vary rapidly with dis- 
tance and angle is not shared by straightforward numerical integration schemes. Based on 
the diffusivity in water of a body of radius 50 nm, we estimate that our basic simulation 
timestep (an average of one attempted virtual move per particle) corresponds to ^ 10 -6 s. 
Our longest simulations exceeded 10 7 timesteps, implying that we probe 'real' timescales in 
excess of a second. 

Our results are summarized by the 'kinetic phase diagram ' of Fig. [J a). This dia gram 
identifies those self- assembled products, whether equilibrated or kinetically trapped, acces- 
sible to dynamical simulation starting from well-mixed initial conditions. We expect by 
extension that such products will be accessible to experiments starting from similar condi- 
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tions. The horizonal and vertical axes of this diagram label the quantities Ri/Ro and e, 
respectively; peanuts with small (resp. large) hydrophilic lobes are found to the left (resp. 
right) of the diagram. We identify regions of compact clusters (squares; see Fig. [jjb)), non- 
spherical micelles (triangles), spherical micelles (circles; see Fig. [jje)) and bilayers (crosses; 
see Fig. [jjc)). Our classification scheme and a summary of the computational resources we 
deployed are discussed in the Appendix. We find these phases to be localized in peanut 
shape space largely in accord with a simple estimate of peanuts' molecular packing param- 
eters [43l [44] (see Appendix) : we have labeled the diagram with lines of packing parameter 
equal to 1 (the regime in which one expects bilayers); 1/2 (the cylindrical micelle regime); 
and 1/3 (the spherical micelle regime). 

The phase most interesting from the materials scientist's perspective is perhaps the bi- 
layer, because extended two-dimensional structures that do not require templating by a 
substrate are attractive candidates for device fabrication. Our observation that amphiphilic 
nanoscale peanuts can form bilayers, and that bilayer formation is localized to a specific 
region of peanut shape space, will facilitate the experimental search for this phase. We also 
observed bilayers to coexist with faceted polyhedra composed of peanuts arranged with a 
high degree of local order (see Fig. [j](d) and Fig. [6](d)). Like micelles, and unlike extended 
clusters, polyhedra present hydrophilic surfaces to solution and are consequently size-limited. 
We do not know why bilayers and polyhedra coexist: they may represent comparable min- 
ima of free energy, or one may embody a particularly stable kinetic trap. This coexistence 
- which resembles simultaneous micellization and phase separation [39] - will be the focus 
of future work, but its observation nonetheless suggests that experimental realizations of 
this system might achieve multipotent self-assembly, wherein components form coexisting, 
ordered phases of strikingly different symmetry. Some protein complexes appear to achieve 
such multipotent assembly by executing conformational changes [45J. 

Of the other phases observed, micelles are found abundantly for many peanut shapes, 
especially when the thermodynamic driving force for association is large and the forma- 
tion of ordered structures like bilayers is hindered by kinetic traps. Micelle polydispersity 
may be controlled by varying Ri/Ro (Fig. [2]), suggesting a route to the synthesis of size- 
controlled nanometer scale assemblies. In addition, cluster formation is possible when the 
hydrophilic lobe is very small. Peanuts whose shapes lie on the left branch of the line 
e = \Ri — Rq\/(Ri + Rq) in Fig. [j](a) are isotropically attractive spheres, and form close- 
packed, crystalline clusters. Peanuts with small hydrophilic lobes also form aggregates whose 
innermost particles make 12 pairwise contacts, and extended crystalline order is possible even 
when the hydrophilic lobe is moderately sized (see Fig. [jjb)). 

While considerations of packing predict the phases we have found, with the exception of 
faceted polyhedra, they give little insight into the complex dynamics of assembly we have 
observed. We focus on the hierarchical dynamics associated with peanut state (i?i/i? , e, e^) 
= (0.9, 0.2, 5/cbT). First, fluctional micelles appear. These can rearrange in a collective 
fashion to form a flattish, metastable proto-bilayer composed of about 7 close-packed peanuts 
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per face, plus attendant edge particles (see Fig. |3](a), left, circled). Such nuclei do not grow 
immediately, but do so only after an additional collective rearrangement that sees an abrupt 
increase in the number of bilayer-like particles comprising their core (see Fig. [3](b,c)). As 
the nucleus exceeds a critical size, bilayer growth proceeds readily; the largest bilayers we 
observed in systems of 1000 peanuts exceeded 800 particles in size (in exploratory simulations 
of 6000 peanuts we observed bilayers comprising in excess of 2000 peanuts) . The time of the 
appearance of the first bilayer-like nucleus exceeding 50 particles in size is long and drawn 
from a broad distribution (from 20 independent simulations of ^ 2.5 x 10 7 timesteps we 
observed the appearance of such structures in 17 cases; of this set, mean nucleation time 
was 1.2 x 10 7 timesteps, with a standard deviation of 9 x 10 6 timesteps). We have not 
explored in detail the dynamics of bilayer formation at other points in peanut shape space, 
but we present one example in Fig. [4| Here a bilayer and a micelle merge before the former 
grows at the expense of the latter. 

Finally, we note that cluster formation by 'peanuts' whose hydrophilic lobes are of van- 
ishing size proceeds via 'two-step' crystallization [46H49] : first, dense amorphous blobs nu- 
cleate from vapor; second, crystalline order nucleates within these blobs (not shown). We 
observed a similar dynamics (not shown) for crystal-forming peanuts (see Fig. [jjb)) possess- 
ing hydrophilic lobes of moderate size. It would make for an interesting theoretical study 
to determine the dynamics of assembly as one proceeds into the regime of increasing hy- 
drophilic lobe size: what is the nature of cluster formation as the crystal structure becomes 
suppressed or thermodynamically disfavored with respect to the dense amorphous phase? 
Clusters in our simulations generally coalesce upon contact, suggesting that in experiment 
such assemblies would not be soluble, and would not be of great practical interest. 

We have demonstrated that dynamical simulation of a simple model of nanoscale am- 
phiphilic peanut-shaped colloids generates phases of clusters, micelles and bilayers in accord 
with expectations based on simple ideas of packing. We have also found certain ordered 
structures and observed complex dynamical mechanisms that cannot be so explained. We 
expect this model and its behavior to be prototypical of a class of amphiphilic structures 
now being synthesized in large quantities, and therefore to help guide the experimental 
search for new, functional nanostructured materials. Future theoretical work involving this 
model will focus on further quantifying the dynamical pathways observed in this study; 
determining the effect upon assembly of liquid-liquid interfaces [50]; and studying the 
design potential of multi-lobed peanut generalizations. 

Note: In a recent study, Miller and Cacciuto explored the self-assembly of spherical 
amphiphilic particles using molecular dynamics simulation [37]. Interestingly, despite the 
differences in interaction range and particle shape between their model and ours, those 
authors also observed the coexistence of bilayers and faceted polyhedra. This implies that 
multipotent assembly of this nature is not dependent upon fine details of particle-particle 
interactions, but can be understood on more general grounds. Furthermore, it may be 
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that such multipotent assembly has already been seen in experiment: one of the referees of 
this paper pointed out that surfactant bilayers can form icosahedral structures in salt-free 
cationic solution [5T] . One possible origin of these structures is the electrostatic force [52j; 
it is also conceivable, based on our results and the results of Ref. [37], that their existence 
might be explained in purely geometrical terms. 
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APPENDIX 



Angular modulation of the hydrophobic interaction. The factors A in Equation ([T]) have the 
following origin (refer to Fig. [5] for geometry). We assume the strength of the hydrophobic 
attraction on the lengthscales we consider to be proportional to the contacting surface area 
between two hydrophobic bodies [40J. We define two surfaces to be 'in contact' if they lie 
no further than 5 nm apart, reflecting our assessment of the likely range of the hydrophobic 
force. If two spheres of radius R w 50 nm are placed in contact, the plane half- angle 
subtended by the limit of the 'contacting' surface area at the center of either sphere is 
ip = arccos(l — 5/ (2-50)) ~ 18.1° (peanuts are azimuthally symmetric about their orientation 
vector, and so we consider only plane angles, as sketched in Fig. [5]). Hence we argue that 
the hydrophobic lobe of peanut % presents its maximum possible surface area to a similar 
body B if the angle between Si and the vector pointing from Q to the center of B is not 
within an angle ifj of the line joining Q to the nearest surface obstruction. 

We next define the angle # max to be the angle between Si and the line joining Q to the 
largest obstruction provided by each hydrophobic lobe's hydrophilic partner. This angle is 
given by geometrical considerations. The angle between Si and the line joining Q to the 
intersection of hydrophobic and hydrophilic lobes is #i n tersec = vr — arccos(^/i? ), where £ = 
5/2 — (R\ — Rl)/(25) (note that £ may be negative). The angle between Si and the line joining 
Ci to the greatest diameter presented by its hydrophilic partner is # g i r th = 7r — arctan(i?i/5). 
We take # max = min(# intersec , g i r th)j unless the largest diameter of the hydrophilic lobe is 
buried in its hydrophobic counterpart (corresponding to Ri < R and t > 5); in this case 

We take # max = ^intersec- 

To construct the function A(9) we argue as follows. Let 9 be the angle between Si and 
the line joining Q to the center of a similar body B. For 9 = the peanut i presents 
its maximum possible hydrophobic surface area to B. Conversely, the greatest obstruction 
of the hydrophobic lobe caused by the hydrophilic lobe occurs when 9 = 7r, i.e. when the 
body B approaches the peanut from its hydrophilic side. This obstruction can be total, if the 
hydrophilic lobe is sufficiently large, or less than total, if the hydrophilic lobe is small and/or 
well-buried within its hydrophobic partner. We take the smallest possible 'signal' presented 
by the hydrophobic surface of peanut i to the body B to be A min = max(0, 1 — (tt — ^ max )/^). 
Finally, we assume for simplicity that the signal interpolates linearly with 9 between its 
largest and smallest values, i.e A{9) oc — 9 if B lies in the penumbra cast by the hydrophilic 
lobe (when 9 exceeds # max — These arguments imply the following piecewise linear form 
for the angular modulation function in Equation ([!]): 



the nanoscale particles we have studied here, identification of the length and angular scales 
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over which the hydrophobic effect operates reveals that multibody forces between particles 
are neither required nor warranted. 

We note that our estimate of the range of the hydrophobic attraction is rough (see 
e.g. [53J). If our assessment of the range of the hydrophobic force is inaccurate, then our 
results will describe nanoparticles whose size is not the 50 nm assumed here. We have 
not performed systematic simulations using different ranges of attraction, but exploratory 
simulations of peanuts whose hydrophobic lobes' radial attractions have a square well form 
of range equal to the lobe radius can form at least one phase not accessible to particles 
possessing the shorter range of attraction studied in the text. This phase consists of 
micelle-like blobs (of many particle diameters in girth) possessing dense, liquidlike interiors. 
These blobs can grow through their mutual coalescence and subsequent restructuring. We 
expect therefore that small nanoparticles (of girth, say, ^ 5 nm) can form phases that 
larger nanoparticles (of size, say, ~ 50 nm) cannot. 

Peanut volume. The volume presented to solvent by two fused spheres of radii R and Ri 
whose centers are 5 apart Rq\ < S < (Ri+Rq)) is given from geometrical considerations 



Here e = 5/(R + i?i), £ = 5/2- (R\ - R 2 )/(25) : E = Rq + R ± and A = R - R v This 
formula may be used to calculate the fraction of the simulation box occupied by peanuts of 
arbitrary geometry. 

Computational details. To obtain Fig. [j](a) we performed simulations for times sufficient to 
observe what we believe to be steady-state behavior at each location in peanut shape space. 
In the regime Ri/Ro > 1.2 (the micelle-forming regime) this is straightforward: for each 
point in the (Ri/Rq^e) plane we performed one simulation at each of = (5, 6, 7, 8)/cbT 
for ^ 5 x 10 6 timesteps (60 CPU hours per simulation). In the cluster- forming regime of 
Ri/Ro < 0.4 we performed 5 independent simulations at each labeled point for each value of 
eb = (4, 4.5, 5, 5.5)/cbT, also for ^ 5 x 10 6 timesteps each (60 CPU hours per simulation). In 
the intermediate regime 0.4 < Ri/R < 1.2 we found some bilayers to form only after very 
long times. We therefore performed 5 independent simulations of ~ 5 x 10 6 timesteps each 
(60 CPU hours per simulation) at each point for each value of = (4, 4.5, 4.75, 5, 5.5)/cbT, 
and in addition performed three further independent simulations of ^ 2.5 x 10 7 timesteps 
(300 CPU hours per simulation) at each point for each of the five values of e^. Along the 
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line Ri/Ro = 1 we performed three additional simulations of ~ 2.5 x 10 7 timesteps for 
each value of 6b = (6,6.5)/cbT. We can of course not rule out the likelihood of nucleation 
and growth of structures on timescales longer than those we probed, nor rule out the 
appearance of other interesting ordered structures at values of that we did not consider. 
To obtain Fig. [2] we performed 10 additional simulations of ~ 5 x 10 6 timesteps at each of 
the four thermodynamic states shown (60 CPU hours per simulation). To obtain nucleation 
statistics for bilayers similar to those seen in Fig. [3]) we performed 20 additional independent 
simulations of ~ 2.5 x 10 7 timesteps (300 CPU hours per simulation) at thermodynamic 
state (i?i/i? ,e,e b ) = (0.9, 0.2, 5k B T). 

Classification of self- assembled structures. We classified the products of peanut self-assembly 
by analyzing the final configuration of each simulation as follows. We define contacting 
peanuts to be those possessing a pairwise energy of interaction of —k^T or less. We define 
an 'aggregate' to be a contiguous set of contacting particles. We identified locations in 
peanut shape space to contain 'clusters' (compact structures with some degree of close 
packing) if any particle possessed 12 contacts. We identified micelles by calculating for 
each aggregate the number Ai = TV" 1 Yli^i ^ ' ( r CM _ Tj). Here N a is the number of 
peanuts comprising the aggregate, Vqm is the center of mass of the aggregate, and V{ is the 
position of aggregate constituent i (all coordinates are corrected for periodic boundaries). 
If any location in shape space contained at least three aggregates of size three or greater 
that scored Ai > 0.9 then we considered that state to contain spherical micelles, and 
marked that state with a circle. Any location with three or more aggregates of size 
three scoring 0.5 < A4 < 0.9 is considered to contain nonspherical micelles, and was 
marked with a triangle. Bilayer order was identified on the basis of the order parameter 
B = 2(N di (N di — l)) _1 ^(ij)(^ ' ^j) 2 -> where the sum runs over all particle pairs in the 
aggregate (JB — | for collections of randomly-oriented peanuts). Configurations possessing 
an aggregate of (25,50,200) constituents or more having a score of B > 0.4 were marked 
with an (orange, red, purple) cross. In Figs. [3] and [i] we show for certain aggregates 
the number of their constituents possessing local bilayer-like order, n^. To compute this 
quantity we formed the dot product of the orientation vector of a given particle and each 
of its neighbors' orientation vectors. We define a 'bilayer-like' particle as one possessing 3 
or more neighbors with which it has orientation-orientation dot product > 0.8, and 1 or 
more neighbors with which it has dot product < —0.8. This order parameter is not useful 
for examining clusters, whose internal particles associate in an orient ationally disordered 
way and can test 'false positive' for local bilayer ordering. However, it does successfully 
distinguish between regions of local order and disorder in the bilayer-forming regime of 
peanut shape space. For most conditions considered it is possible to find amorphous 
structures answering to none of the above descriptions. This is most often true when the 
hydrophobic driving force is very strong, and kinetically frustrated aggregates form. We 
have ignored these structures, focusing on the ordered products we have described. While 
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the order parameters and the threshold numbers we have applied are arbitrary, and the 
extent of each regime varies as the structural classification criteria are varied, we consider 
the trends we have identified to be representative of the nature of the self assembly we 
observed. Some additional structures are shown in Fig. |6| We used VMD [54] to render 
simulation configurations in this paper. 

Packing parameter. Considerations of packing have been used with success to predict the 
phases formed by molecular surfactants [43, 44J. We can define a packing parameter for an 
aggregate as P agg = Kgg/ (^agg^agg), where V agg is aggregate volume, a agg its surface area, and 
L agg a characteristic length. For a cylinder of radius R and length t, for example, Vender = 
7iR 2 £, L cy i inder = R, and a cy i inder = 2ttR£, giving P cy i ind er = 1/2. Analogously, P sph ere = 1/3 
and Pbiiayer = 1- If we now assume such shapes to be composed of amphiphilic components, 
then, loosely, a agg corresponds to the total area presented to solvent by hydrophilic groups, 
while Va gg and L agg are the volumes and characteristic lengths of the hydrophobic units. 
If one assumes that P agg can be estimated from the geometry of a single amphiphile (see, 
however, Ref. [44J), then values of P S i ng i e = P near 1, 1/2 or 1/3 suggest self-assembly 
of bilayers, cylindrical micelles or spherical micelles, respectively. We can make a rough 
estimate of P for an individual peanut. If we assume that hydrophobic and hydrophilic 
lobes touch but do not interpenetrate, then V = ^ttRq, L = 2R and a = nR\. We thus 
have P = |(i? /i?i) 2 . To this approximation, therefore, contours of constant P are vertical 
lines on Fig. [ijb). Furthermore, P = 1 =^ Ri/Rq = \^/3 « 0.82 and P = 1/3 =>► 
RjRo = V2 « 1.41. From Fig. [lfb) we indeed find bilayers and spherical micelles in the 
corresponding regimes. A more refined estimate based on peanut geometry alone may be 
made by accounting for the volume of the hydrophobic lobe buried within the hydrophilic 
lobe. We have 

4 r R o 

V = —ttRq - ix \ dx{R\ - x 2 ), (4) 

L = R + 1, and a = nR\ (valid when the hydrophilic lobe's greatest diameter is not buried 
within the hydrophobic lobe; if this is not true, then a = n(Rl — £ 2 )). Contours of P — 1, 
P — 1/2 and P — 1/3 to this approximation are plotted in Fig. [jja). 
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FIG. 1: Kinetic phase diagram identifying products of self-assembly, (a) Kinetic phase diagram in 
the space of e = 5/ (Rq + R\) and Ri/Ro identifying the products found following dynamic simula- 
tions of peanuts with specified geometries (see text for classification of bilayers, micelles and clus- 
ters). Examples of such products: (b) crystalline cluster at thermodynamic state {Ri/Rq, e, eb) = 
(0.5, 0.5, 4/cbT) (peanuts shown reduced in size); (c) bilayers at (Ri/Rq, e, eb) = (0.9, 0.2, 5/cbT); 
(d) faceted polyhedron at state (i?i/i?o, e, eb) = (0.96, 0.04, 4.75/cbT) (peanuts shown reduced in 
size); and (e) spherical micelles at (Ri/Rq, e, eb) = (2,0.8, 8/cbT) (bottom right shows hydrophilic 
lobes reduced in size). Contours of constant packing parameter (values indicated) are dotted lines; 
the dashed lines indicate the two branches satisfying e = \R\ — Rq\/(R\ + Rq). Below these lines 
one lobe of the 'peanut' is completely buried within the other. 
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FIG. 2: Micelle polydispersity. Number of aggregates of size s, for micelle- forming peanuts 

satisfying e = 1/2 and eb = 8/cbT. Varying Ri/Rq, which is achievable synthetically, controls 
micelle polydispersity. Each line was obtained using 10 independent simulations evolved for ~ 
5 x 10 6 timesteps. 
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FIG. 3: Bilayer self-assembly from one simulation at thermodynamic state (Ri/Rq, e, et>) = 
(0.9, 0.2, 5/cbT). (a) Time-ordered configuration snapshots. Left (t ~ 2.1 x 10 6 ): a micelle es- 
tablishes a bilayer-like configuration, but fails to grow and eventually evaporates. Bilayer-like 
particles (see Appendix) are shown in pink. Center (t « 6.8 x 10 6 ): a similar structure has under- 
gone a collective rearrangement and attained a critical size, and as a result grows readily. Right 
(t « 1.6 x 10 7 ): it is eventually joined by a second bilayer; these grow until they possess about 800 
of the 1000 peanuts in the simulation box. (b) As a function of number of simulation timesteps t 
we plot number of particles n and number of bilayer-like particles rib comprising the structure seen 
in the center panel (which occasionally merges with the second bilayer seen in the right panel). 
Inset: early-time behavior. These quantities are plotted in (c) with t as a parameter. 
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FIG. 4: Bilayer self-assembly from one simulation at thermodynamic state {R\/Rq, e, et>) = 
(0.8, 0.25, 4.75fc B T). (a) Configurations at t = (4,4.5,9) x 10 6 (left to right): a bilayer and a 
micelle merge, with the former eventually consuming the latter. Bilayer-like particles are shown 
in pink. Only aggregates possessing more than 5 constituents are shown, (b) Size n and number 
of bilayer-like particles comprising the most bilayer-like aggregate in the simulation box, as a 
function of t. 
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FIG. 5: Peanut-peanut interaction potential, (a) Peanut geometry with hydrophobic (resp. hy- 
drophilic) lobe colored blue (resp. yellow), together with the radial component of the blue-blue 
interaction potential, (b) Angular modulation function, A(Q), designed to account for the degree 
of contact made by hydrophobic lobes of adjacent peanuts (see Equation [2]). 9 is the plane angle 
between the peanut orientation vector and the vector pointing from the center of the hydrophobic 
lobe (blue) to the center of a similar body, B (shown reduced in size). We plot A{9) for peanuts 
of geometry 1) (i^/i^e) = (0.8,0.85) and 2) {Ri/R^e) = (0.1,0.85). 'Penumbra' and 'umbra' 
identify the regions for which, respectively, partial and total occlusion of the blue lobe is caused 
by the yellow lobe. See text for details. 
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FIG. 6: Additional configurations, (a) Micelles of various morphologies found at thermodynamic 
state (Ri/Rq, e, et>) = (1, 0.8, 5.5/cbT); (b) coexisting bilayers and micelles at state (Ri/Rq, e, eb) = 
(0.8, 0.5, 5/cbT); (c) disordered wormlike micelle at state {R\/Rq, e, eb) = (0.6, 0.8, 5/cbT); and (d) 
coexisting polygon and bilayer at state {R\/Rq, e, eb) = (0.96, 0.04, 4.75/cbT). 



